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Recent studies in ecology and epidemiology indicate that it is important to include spatial 

heterogeneity, synchronization and seasonality in the theoretical models. In this work, spatial het- 

p^ ' erogeneity is introduced via coupled map lattices (CML) and partial differential equations. Stability 

and persistence of some realistic CML are discussed. Chaos control, synchronization and persistence 

are studied for some CML. Some applications in population biology and ecology are given. A simple 

1-^ ' method for finding the sufficient conditions for the existence of periodic solutions for differential 

^r^ , equations with periodic coefficients is given. This will simplify the study of seasonality. 
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I. INTRODUCTION 



Recently it has been argued [|10|,^ that synchronizing local epidemics increases the probability of its global fade 

i-rt out or eradication through careful vaccination program. Similarly synchrony between subpopulations may be an 

PJ important element in species extinction. Therefore studying spatial synchrony, chaos control to achieve synchrony, 

O and persistence in spatial systems is an important task. Also it was shown |10| , p^ that seasonality is an important 

, ^, , factor that should be included. This work is an attempt to model some of these factors simultaneously. 

Chaos is an interesting topic with applications in many fields M,E9l . Most interests are to study the time evolution 

7—i , of systems, while in many realistic systems spatial effects should be included. Hence both cellular automata (CA) 

J> H and spatiotemporal chaos (STC) has arisen [p^p|]. In some cases they generalize differential equations. Typically 

O^ , STC arises in systems with large sizes and in systems consisting of many different components. These two conditions 

Jy^ • are satisfied in many biological and economic systems, hence studying STC in these systems is important. 



1^ The concept of coupled map lattices (CML) was introduced by Kaneko |I4| as a simple model with the essential 

f— ^ i features of STC. It can also be used to approximate partial differential equations. Hence it will be used here to model 
y—( ' some biological and economic systems. In other cases partial differential equations will be used directly to model 
f^ , spatial heterogeneity. 

^^ ■ The paper is organized as follows: In section 2, some applicable CML models are given. A CML for a monopoly 

J^ ' of a company producing different products is studied. Also we have studied some problems in population biology like 
r^ ■ prey-predator, migrating salmon in a river and Schistosomiasis. Chaos control and chaos synchronization in CML are 
discussed in section 3. In section 4, persistence in CML is studied. Section 5 is devoted for developing a method to 
find the sufficient conditions for the existence of periodic solutions for differential equations with periodic coefficients. 
Some conclusions are given in section 6. 



II. SOME APPLICABLE CML 

In this section, some applicable CML are presented, then chaos control and synchronization for CML are studied. 
We begin by a model which is relevant to both economy and biology. Consider a monopoly model where a company 
produces quantity q{t) of a goods, then the profit function is given by 

U{t) ^ q{t)[a ^ bq{t) - c], (1) 

where a, b are some constants and c is the cost per unit of the product. The quantity produced in the next time step 
is expected to be in proportion to the profit i.e. q{t + I) = a H(t), a is a proportionality constant. Rescaling q{t), 
one gets the following dynamical system 



qit + l) = rqit)[^-m] (2) 

where 

q{t) ^ q{t), r = a{a-c). 

a ~ c 

Now assume that the company produces n different products qi, i — 1,2, ...,n (or it has different production sites), 
then the system (2) is replaced by the following CML 

q{i, t + i) = n{i, t) + D [ii{i + 1, t) + n(i - 1, i) - 2n(j, t)] , 

n{i,t)^r{i)q{i,t)[l-q{i,t)], ^^^ 

where the tilde has been dropped for simplicity and D is a, constant representing diffusion of the profits between 
different components of the company. Here D is considered to be small. It is straightforward to find the steady states 
ofEq. (3) 

q{i) ~q, + Dq, , q, ^ 1 - — -, q, = -— !— . (4) 

r[t) r[i) — 1 

To determine the chaotic behavior of the system, local Lyapunov exponents A^, i = 1, 2, ..., n should be calculated. 
It is important to recall that the system exhibits STC if Dl is proportional to the volume of the system, where 

K K K+l 

Dl^K+- i^A., ^A. >Oand ^ A, <0. (5) 

In 1-dimensional systems this is equivalent to the condition that the number of positive Lyapunov exponents is 
proportional to the size of the chain. For D <ti 1, local Lyapunov exponents are given by 



1 ^ 
lim J2ln\{l-2D)r{i)il-2qii,t))\. (6) 

— >oo / ^ — ^ 



t=l 



It is clear that the effect of diffusion is to stabilize the system. This may be interesting economically. 

Consider single species living in patches z = 1, 2, ..., ti with small number diffusing between the neighboring patches, 
then the corresponding CML is 

q*+i = r,ql (1 -ql)+D {qU + ll+i - '^qf) (7) 

The nonzero steady states are given by Eq. (4), and Lyapunov exponents are given by Eq. (6). Again diffusion has 
a stabilizing effect. This confirms that, in general, spatial diversity stabilizes the biological system. 

A more realistic formulation for CML (Eqs. (3), (7)) which takes into account the fact that interactions have finite 
propagation speed is [gl| 

ql^'=nql{l-qD+D{qlzl+qt^l-2qD 

The steady states are not affected by the delay. 

Appling the above analysis to the 2-dimensional Lotka-Volterra predator (y) , prey (a;) CML J20| 

xt+i{i,j) = fJ'XtiiJ) [1 - xt{i,j) - yt{i,j)] + 

Di [xtii + lj) +xtii- l,j) +xt{i,j + l)+xt{i,j - 1) -4xt(«,j)], 

yt+i{hj) = Pxt{i,j)yt{iJ) + 

D2 [yt{i + 1, j) + ytii - i,i) + yt{i,j + i) + yt{i,j - i) - 4yt(i, j)] , 

where fj,, /?, Di and Di are some positive constants. The coexistence steady state (a; 7^ 0, y 7^ 0) are 

, , 1 , . 11 

The largest local Lyapunov exponent is 



(8) 



X{i,j) =limT^ooTELol^l/t(*'i)l 



ft{i,j) ^ \ \yt{i,i) + y/U?{i,j) - APfiXtii,j)yt{i,j) 
Ut{i,j) - A^ (1 - 2xt{t,j)) + (3xt{i,j) - 4(Di + D2) . 



(9) 



Proposition 1: If 1 < r{i) < 4 Vi = 1, 2, ...n, then the interval [0, 1] is a local trapping region for the systems (3) 
and (7). 



Proof. If Iqjl < 1, then 



k?+i < (1 - 2D) HI +D w+' + nj-M < 



nil -2D) ^ D{n^i + n+i) 



< 1. 



Also using the results of [|], then it is direct to prove the following proposition. 

Proposition 2: If D is sufficient small and r{i) > 4, Vi = l,2,...n, then the symbolic dynamics of the systems 
(3) and (7) is approximated by the direct product of the symbolic dynamic of the local subsystems. 

An interesting application of CML to a realistic ecological problem is the problem of migrating Salmon in a river 
|L8| . Typically the predators occupy fixed positions and their number is approximately constant, hence the system 
can be modelled by the following equations 



^t+i 






(10) 



where xl is the number of salmon at site i and time t, y* is the number of predator at site i (assumed to be time 
independent), and 7 is a positive constant. The number of predated salmon in the i-th patch is 72/V(l + ^^t); hence 
it decreases as the patch size (a;J) increases. This agrees with observations. The steady state solution of (10) is 



1 ' r / 

^' = ^n V^a;'=-i + (l + 7y'=)2-(l+7/) 



fe=i 



(11) 



ej. Then it is straightforward to obtain 



To study the stability of Eq. (11), one imposes a small perturbation xl — 
that el — X]fc=i [~7y'^/(l + 2;'^)] • Since, typically, x'' » y'"', then we have: 

Proposition 3: For all finite values of n, the steady states (11) are stable. 

Now CML for an epidemic model is constructed. Consider a certain population, according to the health of each 
individual, the population is classified to susceptible (S), infected (I) or recovered (removed) (R). Assuming that only 
infectives (and infected) diffuse, then a CML for the susccptible-infected-removed-susceptible model [Q (SIRS) can 
be written as 



4-f3sm+aRi, 
J P „ ^, P 






'il,^n+psiii~^ii+D(ir 



lIl-aRl 
211), 



P-^ 



(12) 



where a, P and 7 are positive constants. 

For the realistic case of rabies, the CML model can be written in the following form: 



J J 



f3sm, li^,^ll+Pslll-^ll+D{ll+'+ll 



2I{). 



(13) 



The largest local Lyapunov exponent is given by 



lim 



1 ^ 



t=i 



a + VQ 



Ab 



where 



a == 2 - 7 - 21) - /3/f + Ps{, 6 = 1 - 2D - 7 - /3//(l - 2L> - 7) + /?s{- 

Now the CA'EL of schistosomiasis is discussed. Schistosomiasis is one of the major communicable diseases of pubhc 
heahh. It has socio-economic importance in many developing countries. Moreover, it has been observed that infection 
is concentrated in certain regions (centres), hence the dynamics of the disease depends on the interaction between 
infected persons at these regions. This can be represented as CML. We begin by studying the dynamics at one 
centre, then later introduce the coupling between them. Field observations have indicated that there are two periods 
(May- June and October-November), where infection of humans by cercaria is highly probable. This makes difference 
equations more suitable to model schistosomiasis with one time step approximately equal six months. The basic 
assumptions are: 

(i) Human population grows exponentially. 

(ii) Total snail population is approximately constant. 

(iii) Infected snails do not reproduce. 

Thus the local equations can be written as 

Ht+i = aHt - (iPu Pt+i = 7Pt{l - ^) + XCtHt, Ct+i ^ tiPt ~ vCu (14) 

a 

where Ht is the human population at time t, Pt is the parasite (adult worm) population at time t, Ct is the cercaria 
population at time f, and a, /?, 7, A, /i, v^ d are positive constants. To keep in touch with observations, we will not 
scale away any of them. The steady states of this system are 



P = or P 



1/7-1 



(15) 



AM/3/[7(i.+l)(a-l)]-l/d' 



where a>l, 0<7<1. These conditions agree with observations, since populations in developing countries are 
increasing. The condition < 7 < 1 means that parasite population increases only through infections by new cercaria. 
The stability analysis depends on the eigenvalues of the Jacobian matrix A 



A 



a -(3 

\C 7(1 - f ) \H 

/i —v 



The zero solution is unstable, while the nonzero solution is stable if the following conditions are satisfied: 

(^+l)(a-l)" 

7(l + a)(l-^)+(,;^§f^>0, (16) 



«7K1 - ^ 



2P_ _|_ apXt^P _ upXfiP 



< 1. 



The first inequality is the familiar threshold condition. 

There are two ways to resist schistosomiasis namely treatment and snail fighting. Both decrease the parameter A 
in Eq. (14). If the treatment is perfect, then A^O, P = C = 0, a ^ 1 and H = constant. However it is known that 
no treatment is perfect, hence the maximum allowed value for A is 

^^ 7(z.+ l)(a-l) 

Allowing the infected persons to diffuse, one gets the following CML 

Hl^, = aHl - f3Pl + D{Hl-' + Hl+' - 2HI), 

Pl^^^^Pl{l-^) + \ClHl (17) 

Cl+,=iiPl-uCl. 

In the limit 13 <C 1, the steady states are given by Eq. (15). From the above discussions, we conclude the following: 



(i) Treatment of most infected persons should be done simultaneously to prevent reinfection through the diffusion 
between sites. This agrees with the idea of synchronization mentioned before. 

(ii) Snail fighting should be done as efficient as possible, since they are crucial in the disease spread (7 < 1). 

(iii) More than one drug (with independent mode of action) should be used alternatively to avoid drug resistance. 

An alternative model to the ones used so far to study spatial effects in epidemics is given by p5| 



dh 
dt 






where Si(Ii) is the fraction of susceptible (infected and infective) individuals in the i-th patch, i — 1, 2, ..., n, and 7, fi 
are constants, /3y are a measure that an infective from the j-th patch will infect a susceptible from the i-th patch. 
The total population in each patch is assumed to be constant. It is realistic to assume that 

Pit > Pij , i 7^ j- 

Expanding Si = s? + sj, h = l'^ + 1} and linearizing in /3y, i ^ j, s\ and 1} , one gets the following equilibrium 
solutions 



A, 



f3i^{^J' + l) 



Y.j^iPi]Ij 



^i + (3^^I°-^lS°,P^^/{^l + J) 



, I 



-fJ- 



M + 7 



It is clear that if the patches are different (i.e. Pa changes significantly with i), then the equilibrium solutions are 
asynchronous. Using the assumption that Pa 3> Pij, i ^ j, the approximated solution becomes 

Si = I, Jj = Vi = l,2,...,n. 

This solution is asymptotically stable if the following determinants are positive 



ai 1 
03 a2 



ai 


1 
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02 


Ol 


05 


04 


as 



.A„ 



ai 1 

as 02 ai 

as 04 as 





where 



at^ Pa- ^1- 7. 



In this work, terms were kept to the second order in Pij, i ^ j. 



III. CHAOS CONTROL AND CHAOS SYNCHRONIZATION IN CML 



In the previous sections, it has been shown that chaos exists in some realistic CML. In some cases chaos is not 
an acceptable phenomena, hence chaos control is required. Chaos control in ordinary systems has been discussed in 
several references e.g. [EgUI^JI]. To control chaos in CML, we propose to increase the coupling constant D. Using 
numerical simulations we found that in a logistic CML (Eq. (3)) with n = 6, such that only one of them with 



r{i) — 2 (the other five oscillators had r{j) — 4) was suSicient to control the chaos of the whole CML provided that 
0.5 > L» > 0.493. 

In some cases it is not required to control chaos, but to synchronize it [f7Uq|. In the work of Codreanu and Savici 
, they synchronized chaos of two identical logistic oscillators using a third one. Their work has a drawback that 
there is no bound on the used control. In many cases this is not realistic. Here their work is generalized to CML and 
bounded control is used to synchronize chaos. Consider the logistic CML (Eq. (3)), with r{i) = 3.7 + 0.3*rnd, where 
rnd is a uniformly distributed random number between zero and one. Thus all oscillators are chaotic but different. 
The symmetric case can be studied similarly. Assume that it is required to synchronize the CML to the a certain 
oscillator (without loss of generality call it the l-th oscillator) . Then to each oscillator i ^ I apply the following control 

U if |u| < Umax, 

U*^^ = { "max if W > Umax, (18) 

"^niax 11 W, -C "~timax5 

where 

u = -r(z)g* + rmqlf + r{l)q\ - r{mY + \{4 - <il), 

and Umax is the maximum allowed control. We found numerically that if n = 10, then Umax — 0.09 is sufficient to 
synchronize the CML. Furthermore the coupling constant D was irrelevant to the value of Umax, which is contrary to 
the case of chaos control. 

Proposition 4: There is bounded Umax, such that the control (18) will synchronize the CML (3). 

Proof. Since g* in Eq. (3) is bounded, then 3 (/max < oo, such that g* < gmax- Choosing Umax = 9max, 
then 

u*+l = -r(z)<z* + riiMf + r{l)ql r{l){qlf + \{4 - q\) if \u'+'\ < Umax. 

Thus the controlled deviation 

A*+^ = 9* - 9? + "I 
satisfies 

A*+i = Ia* i.e. A* ^ 0. 



In practice, simulations have shown that the required Umax <C ^max- The previous study is applicable to any CML. 

IV. PERSISTENCE IN SOME SPATIALLY INHOMOGENEOUS SYSTEMS 

Typically a 1-dimension CML is given by 

0*+i = (1 - D)el + |[0*+i + 0*_j + /(0*), (19) 

or 

0f 1 = (1 - D)fi9l ) + f [/(0*_i) + /(0*+i)], (20) 

where t — 1,2, ... , j — 1,2, . . . ,n. 

Definition 1: A dynamical system is persistent if Va;^(0) G int(S'), then limf^oo inf .^^(t) > Vi = 1, 2, . . . , n. A 
CML is persistent if all its components are persistent. 



An open question is what is the effect of diffusion on persistence? 

Proposition 5: (i) If 3e > such that 'iO e [0,e], then f{e) > 0, and {1 - D + f'{9)) > 0, then the CML 

(19) is persistent. 

(ii) If f{9) > 0, ¥6* > and 3e > 0, such that 6 e [0, e] => f'{e) > t] > 0, and 77(1 - D) > 1, then the CML 

(20) is persistent. 

Proof, (i) V0*- G (0, e], using mean value theorem, then 30 G [0, ^J], such that 

0*+i - (1 - D)9l + ^[9l+, + 9l_^ + f{(:])9] > (1 - D)9] + ^[9]+, + 9]_^]. 



Hence 



9f'><t>f'^{l-D)9] + ^\9]^,+9]_,]. 



2 
The right hand side is a circulant matrix whose eigenvalues are known, and its maximum eigenvalue is unity. Thus 

<^*+i = 0° = 00 ^ 9]+' > 9° > 0. 

(ii) For CML (20), the conditions of the proposition imply 9*-^^ > (1 — D)f{9*j). If 0* G (0, e], then using mean value 
theorem one gets (so long as 6^ G (0, e] Vs < t, otherwise 9^^ > e which implies persistence), i.e 

9f^ > (1 - D)r]9] > [(1 - -0)77]*^" > 0" > 0. 



To study persistence in partial differential equations, we use the following results |^] with slight modifications: 

Theorem 1: (i) Suppose that /(T^,m) is Lipschitz in "f* and continuously differentiable in u with 

df _, 

Te- < for u > 0, /(T", u) <Oiiu>L 

ou 

for some constant Z, and f{~r',Q) > at some point in the domain fJ. Then the following problem 

— = DV^u + uf{T',u) inf7a;(0,oo), (21) 

with Dirichlet or Neumann boundary conditions has a unique positive steady state u"* which is a global attractor 
for nontrivial non-negative solutions (hence the system (21) is persistent), if the following problem has a positive 
eigenvalue a 

au = D\7'^u + uf{-r',0) in f], (22) 

with the same boundary conditions as (21). 

(ii) Suppose that /i(^,ui,U2) and /2(^, 1*1,1*2) sue C^, bounded that /i(T^,ui,0) satisfies the conditions of part 
(i), with positive a, and /2(~?^, 0, U2) < for U2 > 0. Let cto be the largest eigenvalue of the system, then 

aQUi= DiS/^ui+uifi{'r,0,0) inf7, (23) 

with Dirichlet or Neumann boundary conditions. Also let (72 be the largest eigenvalue of the system 

a2U2 = D2\7^U2 +U2f2{^,ur,0) inf], (24) 

where wf^ is given by fi{l^, w^*', 0) — 0, then the following system is persistent if both ao and (T2 are positive 



^ = Dl^^Ul + Ulfli^ ,Ul,U2) 
^ = D2V^U2 + U2f2{^,Ul,U2) 



a% : n\72r : r. )4'r':'(' inf^:.(o,oo). (25) 



Now theorem (1) is applied to some 2-dimensional systems. The domain fl is chosen to be the rectangle x S 
[0, Li], y g [0,^2]. An example for two species problems is host-parasitoid, where ui (1*2) is host (parasitoid), 
respectively and /i and /2 in Eq. (25) are given by 



and persistence conditions become 



/l = 


1-aui ui 


" I + U2 ^' ■" I + U2 




^>l + i?2(fj + f^), 




b-l>Di{^ + ^). 



1, (26) 



(27) 

Another example is the predator-prey or epidemic model where ui (1x2) is prey (predator) or susceptible (infected 
and infectious) and /i and /2 in Eq. (25) are 

fl = l-Ui-U2, f2=Ui-J, (28) 

and persistence conditions are 

22 22 

l>^i(^ + ^)andl>7 + i^2(^ + ^). (29) 

For the case of rabies one sets Di = 0, since only infected foxes spread, while susceptible ones do not. This indicates 
that a method to fight infection spread is to put susceptibles in small groups (to decrease the domain). 
Now the Krill- Whale model [Q will be studied including diffusion. The functions /i and /2 in Eq. (25) are 

Jl ^ rUi{l - Ui) - U1U2, /2 = S-U2 , 

Ul 

where ui (U2) represents krill (whale), respectively and r, s are positive constants. The uniform coexistence steady 
states are mi = 1 — l/r, M2 = 1, r > 1. Its stability requires s < 1 — 1/r, 1 < r < 2. Sufficient conditions for persistence 
can be derived (using theorem (1)) in the following form 

22 22 

r>D,i^ + ^) + l, s>D2{^ + ^) (30) 

^i ^2 ^i ^2 

Furthermore using Dulac-Bendixon lemma, one can see that this system has no cycles. 

There is a weakness in the definition of persistence in that it does not guarantee that the population is larger than 
a prescribed value. To find the minimum of the population is part of the practical persistence [0;p|. Studying this 
concept for the krill-whale model, one can derive the following limits for the populations 

1 > Ml > 1 , U2 — 1. 

r 

These values agreed with the numerical simulations. 

One may argue that advection should be included in the above diffusion model. So an interesting question arises 
as to what is the effect of advection on the persistence conditions of theorem (1)? Consider the following diffusion- 
advection system 

^DV'^u + uf{-t,u)-v— inf7x(0,c»). (31) 

at ox 

If the system (31) admits a solution of the following form 

/ \ v^ \ f/ ■ nnx mix. ,„„, 

M(x,i) = > e " (a„sm— ho„cos— — ), (32) 

^ Ij ij 

and linearizing around m = 0, one gets that the advection term in Eq. (31) affects only the imaginary part of the 
eigenvalue a in Eq. (22). Hence we conclude the following proposition: 

Proposition 6: If the system (31) admits a solution of the form (32), then the persistence conditions of theorem (1) 
are valid. 



V. A METHOD TO STUDY DIFFERENTIAL EQUATIONS WITH PERIODIC COEFFICIENTS 

We begin by proposing the following correspondence between differential equations with periodic coefhcients and 
some difference equations. Consider the following system 



dt 
Approximating the left hand side using 



-TT^Yl P'3 (*)^J' ' P'l (*) "" P^3 (t + T). (33) 



dxi Xi{t + T) — Xi{t) 



dt T 

thus one gets 

x^{t + T)^x,{t)+T^p,,{t)x,, t = 0,T,2T,.... (34) 

j 

Since the discrete system (34) has time step T, then Pij{t) will be constant in its evolution, thus the system (34) is 
autonomous. 

Recalling that the stability conditions for discrete systems are more stringent than the corresponding continuous 
one, then stability conditions for the system (34) will be sufficient conditions for those of (33). Thus we have 

Proposition 7: If all the eigenvalues A of the matrix [5ij +Tpij{t)] satisfy |A| < 1 Vt e [0, T], then the periodic 
solution of Eq. (33) is stable. 

As an example consider the following equation 

'Ll+p{t)x^Q, p{t + T)=p{t). (35) 

Approximating it by 

x(t + T) + x(t-r) + a;(t)(-2 + TVi)) = 0, i == T, 2r, 3r, . . . . (36) 

Notice that in Eq. (36), p{t) is constant since t increases by T, and p{t) = p{t + T). Thus one gets that Eq. (36) is 
stable if 

0<T2p(t) <4Vie [0,r], (37) 

which yields the known Lyapunov theorem for the sufficient conditions for the stability of the periodic motion of Eq. 
(36) 

p{t) > 0, < T /" p{s)ds < 4. (38) 

Applying these results to Lame's equation 

f^-{h + 2k\nH)x = 0, 
dt-^ 

one gets the following sufficient conditions for the stability of its periodic motion 

nH) < 4, yt 



< -^K'^{h + 2PsnH) < 4, Vt e [0, 2K], 

^ ^ j^/2 d4> (39) 



Applying them to Hill's equation 



--TT + (a - 2bcos2TTt)x = 0, 



one gets that the sufhcient conditions for the stability of its periodic motion are 



0< a-26cos27ri < 4Vt G [0,1]. (40) 

The previous method is also useful for the case of periodically forced Lotka-Volterra predator prey models of the 
following form 

dp dv \ / ^ \ 

—- — —ap + ppv, —- = [r + esinzTrtjvn ) — ^pv. (41) 

dt dt A 

Following the above procedure, one gets the following stability conditions: 

r + e sin 27ri , , r , , ^ 

0<a — < 4 Vie [0,1]. (42) 

pK 

This method is expected to simplify the study of seasonality. 

VI. CONCLUSIONS 

In conclusion, ecologists have stated that it is important to include several factors in some biological and epidemic 
models e.g. spatial heterogeneity, synchronization and seasonality. Since a crucial question is whether the epidemic 
(or the population) will persist, we studied the persistence of some spatially heterogenous systems (using CML or 
partial differential equations). Synchrony in CML and CA is studied. Since chaos implies asynchrony, we studied 
bounded chaos control on CA and CML. A simple method to study seasonally in forced systems (systems with periodic 
coefficients) is introduced. 
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